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We study pre-equilibrium giant dipole resonance excitation and fusion in the neutron-rich system 
"2Sn-F*«Ca at energies near the Coulomb barrier, and we compare photon yields and total fusion 
cross sections to those of the stable system ^■^''Sn-|-'"'Ca. The dynamic microscopic calculations are 
carried out on a three-dimensional lattice using both the Time-Dependent Hartree-Fock method 
and the Density Constrained TDHF method. We demonstrate that the peak of the GDR excitation 
spectrum occurs at a substantially lower energy than expected for an equilibrated system, thus 
reflecting the very large prolate elongation of the dinuclear complex during the early stages of 
fusion. Our theoretical fusion cross-sections for both systems agree reasonably well with recent data 
measured at HRIBF. 

PACS numbers: 21.60.-n,21.60. Jz 



I. INTRODUCTION 

Heavy-ion reactions at radioactive ion beam (RIB) fa- 
cilities enable us to synthesize new exotic neutron-rich 
nuclei far away from the valley of stability and to study 
their physical properties. At the HRIBF facility a se- 
ries of experiments has been carried out with radioactive 
^■^^Sn beams on various targets ^rt4i]. 

One of the major open questions in fusion reactions of 
exotic neutron-rich nuclei is the dependence of the fusion 
cross section (and related observables) on the neutron ex- 
cess, or equivalently on the total isospin quantum number 
Tz = {Z — N)/2. To reveal possible systematic trends re- 
quires both theoretical and experimental studies with a 
wide variety of projectile and target combinations which 
are expected to become available at current and future 
RIB facilities. To be able to pin down the isospin de- 
pendence in a fully microscopic theory, it is desirable to 
choose collision partners which are as simple as possible: 
projectile and target nuclei should be spherical in their 
ground state, and the compound nucleus should have a 
high fission barrier so that the fission component can be 
ignored (at least at lower beam energies). Ideal reaction 
partners satisfying these criteria include the neutron-rich 
system ^^agn-f^^Ca and the stable system i24Sn-H40cg^ 
which we investigate in this paper. 

If the reaction partners have appreciably different N/Z 
ratios, the proton and neutron centers of mass of the re- 
acting system may not coincide during the early stages 
of the collision. In this case, collective oscillations of 
protons against neutrons on the way to fusion might oc- 
cur, the so-called pre-equilibrium giant dipole resonance 
(GDR) In fusion reactions, the shape of the pre- 

equilibrium di-nuclear complex exhibits a very large pro- 
late deformation as compared to the shape of the equi- 
librated compound nucleus (see e.g. Ref. [l^l. Figures 
1 and 2). Consequently, the 7-rays emitted during the 
early phase of fusion contain information about the shape 



of the system and provide a cooling effect along the fusion 
path as well as insight into the charge equilibration. 

The time-dependent Hartree-Fock (TDHF) theory pro- 
vides a useful foundation for a fully microscopic many- 
body theory of large amplitude collective motion [ill, El] 
including nuclear molecular resonances, deep-inelastic 
scattering and heavy-ion fusion, and collective excita- 
tions. The earliest TDHF calculations carried out in 
the 1970's and 1980's were severely limited by computer 
memory and speed. Hence, these calculations had to 
make restrictive assumptions such as axial symmetry and 
the omission of spin-orbit coupling. During the past 
several years, with a steady increase in computational 
power, it has become feasible to perform TDHF calcula- 
tions on a three-dimensional (3D) Cartesian grid with no 
symmetry restrictions whatsoever and with much more 
accurate numerical methods [H-lig. At the same time 
the quality of effective interactions has also been substan- 
tially improved [17h21j . Our TDHF code utilizes the full 
Skyrme interaction, including all of the time-odd terms 
in the mean-field Hamiltonian [l^ . 

During the past several years, we have developed the 
density-constrained TDHF method (DC-TDHF) for cal- 
culating heavy- ion potentials [l^l, and we have applied 
this method to calculate fusion and capture cross sec- 
tions above and below the barrier. So far, we have 
studied the systems "2Sn-h64Ni_j23, 64Ni+64]xjj j^l, 
i6Q_^208pi3 [2|^ I32,i24gjj_^962i. [ijfand we have stud- 
ied the entrance channel dynamics of hot and cold fusion 
reactions leading to superheavy element Z = 112 [27| . 
In all cases, we have found good agreement between the 
measured fusion cross sections and the DC-TDHF re- 
sults. This is rather remarkable given the fact that the 
only input in DC-TDHF is the Skyrme effective N-N in- 
teraction, and there are no adjustable parameters. 

This paper is organized as follows: in Section II we 
summarize the theoretical formalism for pre-equilibrium 
GDR excitation in heavy-ion fusion reactions with the 
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TDHF method, and we compare the resulting dipole ra- 
diation yields for the neutron-rich system ^"^^Sn-f^^Ca to 
the photon spectra of the stable system ^^^Sn-|-^°Ca. In 
Section III, we study total fusion cross sections for both 
systems. First, we carry out unrestricted TDHF runs at 
energies above the barrier (no tunneling). Secondly, we 
calculate the heavy- ion interaction potential V{R) and 
the coordinate-dependent mass parameter M{R) utiliz- 
ing the DC-TDHF method. By numerical integration of 
the Schrodinger equation for the relative coordinate R 
with the Incoming Wave Boundary Condition (IWBC) 
method, we calculate total fusion cross sections both be- 
low and above the barrier, and we compare our theo- 
retical results to recent data measured at HRIBF. Our 
conclusions are presented in Section IV. 



II. TDHF STUDY OF PRE-EQUILIBRIUM GDR 
EXCITATION 



In the present TDHF calculations we use the Skyrme 
SLy4 interaction for the nucleons including all of the 
time-odd terms in the mean- field Hamiltonian (l3| . The 
numerical calculations are carried out on a 3D Cartesian 
lattice. For the two reactions studied here, the lattice 
spans 50 fm along the collision axis and 30 — 42 fm in 
the other two directions, depending on the impact pa- 
rameter. Derivative operators on the lattice are repre- 
sented by the Basis-Spline collocation method. One of 
the major advantages of this method is that we may use 
a relatively large grid spacing of 1.0 fm and nevertheless 
achieve high numerical accuracy. First we generate very 
accurate static HF wave functions for the two nuclei on 
the 3D grid. The static HF equations are solved with 
the damped gradient iteration method. The initial sep- 
aration of the two nuclei is 22 fm for central collisions. 
In the second step, we apply a boost operator to the 
single-particle wave functions. The time-propagation is 
carried out using a Taylor series expansion (up to orders 
10 — 12) of the unitary mean-field propagator, with a time 
step M = 0.4 fm/c. 

In TDHF, the dissipation of the relative kinetic energy 
into internal excitations is due to the collisions of the 
nucleons with the "walls" of the self-consistent mean- 
field potential [2^ [2^ . Because the randomization of the 
single-particle motion occurs through repeated exchange 
of nucleons from one nucleus into the other, the equili- 
bration of excitations is rather slow, and it is sensitive to 
the details of the shape evolution of the composite sys- 
tem. This is in contrast to most classical pictures of nu- 
clear fusion which generally assume near instantaneous, 
isotropic equilibration. 

In Fig. [T] we show the time evolution of isoscalar mul- 
tipole moments for a head-on collision of Sn on Ca at 
Ecm = 130 MeV. Compared are the two isotopic com- 
binations i32Sn-h48Ca and ^'^'^^n+^^C&. The isoscalar 
quadrupole moments Q2Q and the associated quadrupole 
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FIG. 1. (Color online) Time evolution of isoscalar multi- 
pole moments for a head-on collision of Sn on Ca at Ec-m = 
130 MeV. Two isotopic combinations are considered as indi- 
cated. 
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show a remarkable stability in the pre-compound state, 
i.e. after 1000 fm/c. However, this extraordinarily long 
life-time may be partially due to the frozen occupations 
(no self-consistent pairing) in the TDHF calculations and 
partially due to other missing correlations. The center- 
of-mass coordinate Xc.m. moves less than 0.01 fm. The 
isoscalar r.m.s. radius shows clearly a small, but steady 
growth. This is most probably caused by nucleon emis- 
sion in connection with reflecting boundary conditions. 
Nucleons which leave the compound system are spread 
over all box volume and deliver unnaturally larger con- 
tributions to the radius. A computation of the radii in 
an analyzing volume confined to the vicinity of the com- 
pound system would probably show the opposite, slightly 
shrinking radii. 

The excitation of the pre-compound giant dipole mode, 
which is likely when the ions have significantly differ- 
ent N/Z ratio, is a refiection of dynamical equilibra- 
tion. Dipole excitation results in the emission of 7 ra- 
diation and possibly pre-equilibrium particle emission 
which in turn provides an important cooling mechanism 
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FIG. 2. (Color online) Time evolution of isovector multi- 
pole moments for a head-on collision of Sn on Ca at Eon = 
130 MeV. Two isotopic combinations are considered as indi- 
cated. 



during the early stages of the fusion reaction. Sev- 
eral transport simulations like semiclassical Boltzmann- 
Nordheim-Vlas ov [ 5, 30] as well as microscopic calcula- 
tions [U, |3ll, nil have confirmed this scenario. 

To calculate the electric dipole radiation emitted by 
the time-dependent charge distribution, we consider the 
dipole operator 



, , NZ 
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which represents the expectation value of the x- 
component of the dipole operator d^/e taken with the 
time-dependent TDHF Slater determinant |$(i) >. 

In Fig.[2]we show the time evolution of the correspond- 
ing isovector quadrupole moment (top panel) and the 
isovector dipole moment D{t) (bottom panel) for the two 
systems. It is interesting to compare the isovector dipole 
moments D{t) for i^sSn-h^^Ca and for ^"^^Sn+^^Ca.. The 
initial amplitude seems to be nearly the same in both 
cases which may question the argument that the total 
isospin Tz is responsible for the strong dipole excita- 
tions. We have also computed a collision of the sym- 
metric system ^®Sr-t-*^Sr which ends in a comparable 
compound state. There is no isovector dipole excitation 
in this case as expected. There is another interesting 
feature in that the dipole amplitude remains rather sta- 
ble for -'^^^Sn-l-'^^Ca while it shows signs of damping for 
i24Sn+40Ca. 

Following the bremsstrahlung approach developed by 
Baran et al. [sO] we define the dipole acceleration 
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and introduce its Fourier transform 
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Alternatively, for nearly harmonic vibrations one can use 



the expression \D"{uj)\'' 



\D(uj)\ , with 
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By adding and subtracting the position vectors of the 
neutrons we can write this expression in the form 
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(3) 

We now impose the condition that the center-of-mass is 
fixed at the origin of the coordinate system at all times. 
As a result, the last term in the round brackets vanishes 
and we obtain the dipole operator in the form 
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Let us now consider a central collision in x-direction and 



D{uj) = / D(t)e'"' sin^ 



t-u, 



dt 



(8) 



The time filtering sin is used to smooth out peaks com- 
ing from finite integration time. The "power spectrum" 
of the electric dipole radiation is given by ^] 



dP 

dK 



2a 



c 



(9) 



where a — /{he) k, 1/137 denotes the fine structure 
constant. 

In Fig. 13] we show results of TDHF calculations of 
the photon yield given by Eq. [Sj For our calculations 
the use of D"{t) or the alternate expression with D{t) 
give essentially the same result. Compared are the 
neutron-rich system ^■^^Sn-|-'*^Ca and the stable system 
^^^Sn-|-^°Ca at the same center-of-mass energy. In per- 
forming the full Fourier transform we have cut out the 
approach phase of the two ions corresponding to times 
i < imin = 800 fm/c. The upper limit of the time in- 
tegration is traax = 7,200 fm/c which corresponds to 
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FIG. 3. (Color online) Power spectrum of pre-equilibrium 
dipole radiation following a central collision at Ecm. ~ 
130 MeV of i^'^Sn-f^Sca (left panel) and ^'^'^Sn+'^^Ca (right 
panel). Results are shown for spectral analysis over the whole 
simulation time (lower panels, red) and over three shorter in- 
tervals as indicated (upper panels). The radiation power spec- 
trum of the (deformed) ground state of the compound system 
is shown by the dashed curve (lower panels). 



2.16 X 10^^° s. The lower part (red curve) of the Fig. [3] 
shows the spectrum obtained for the whole time inter- 
val 800 — 7,200 fm/c. We observe that the isovector 
dipole radiation spectrum shows a peak near 10 MeV 
for both isotopic combinations. Furthermore, we have 
also divided the total time span into three intervals and 
evaluated the power spectra separately for these (see up- 
per panels in Fig. [3l This indicates the time evolution 
of the power spectra. The case of -'^'^^Sn-l-^^Ca shows a 
steady upshift of the mean peak and decrease of overall 
strength. The latter is an obvious effect of damping. The 
upshift indicates a decreasing system size. This complies 
with the expectation that the true system radii (evalu- 
ated in an analyzing volume) should shrink. The sys- 
tem -'^^^Sn-f^'^Ca behaves differently. The spectra stay 
stable during the first two time intervals and shrink dra- 
matically in the third time interval which complies with 
the fact that the dipole amplitude is visibly smaller at 
later times (see D{t) panel in Fig. [2]). In light of the en- 
trance channel N/Z ratios for the two systems, 1.48 and 
1 for ^^^Sn-)-^°Ca stable system and 1.64 and 1.4 for the 
"'^^^Sn-l-^^Ca neutron-rich system, one can conclude that 
the system with a larger initial N/Z differential, namely 
the stable system, should have the stronger initial pre- 
equilibrium GDR excitation. As can be seen from Fig. [3] 
this is indeed the case. On the other hand the presence of 
extra neutrons for the neutron-rich system seem to pro- 
long the life of the GDR suggesting a longer equilibration 
time for this case. 

We have also evaluated the dipole radiation spectra 
originating from the compound systems in their ground 
states. Both of these have a quadrupole deformation 



f^20 ~ 0.3 which is significantly smaller than the defor- 
mation of the pre-compound state, /32o ~ 0.6. Figure [3] 
(lower panels, dashed black curves) shows also the power 
spectra produced by the compound systems. Because of 
the smaller deformation the spectra peak at a higher en- 
ergy, and this peak is well separated from the peak of 
pre-compound stage. 

In order to estimate the energy content in the 
isovector-dipole oscillations, we have computed the 
dipole polarizability ao and find for the pre-compound 
system Z/N = 70/110 the value ao = 30 fm^. The 
corresponding excitation energy is given by Edip — 
e^D^/ {2a]j). This energy content changes based on what 
stage of the reaction is considered. During the initial con- 
tact phase, where the dipole excitation may be the dom- 
inant form of collective excitation the amplitude D (see 
Figure[2]) is as large as 5—6 fm, which yields 0.6—0.8 MeV 
excitation whereas at later stages the dipole amplitude 
is in the range of 2 — 3 fm resulting in approximately 
0.1 — 0.2 MeV excitation. Our calculations of excitation 
energy based on the DC-TDHF method (see Ref. Q) 
for this cm. energy indicate that the amount of intrinsic 
excitation energy starts to rise from zero as the two nu- 
clei make contact and reach the value of 40 MeV as the 
initial composite system is formed. Thus the fraction of 
the excitation energy carried by the pre-equilibrium GDR 
can be a large part of the small excitation energy near 
barrier (by as much as 40% at the barrier peak), while it 
becomes less than a percent in the interior region. 



III. FUSION 
A. TDHF calculations at energies above the barrier 

Depending upon beam energy and impact parameter, 
the dominant reaction channels are deep-inelastic and fu- 
sion reactions [l^l- In general, central collisions and col- 
lisions with relatively small impact parameter result in 
fusion, while larger impact parameters give rise to deep- 
inelastic reactions. At i?c.m. energies above the potential 
barrier it is possible to carry out standard unrestricted 
TDHF runs (no tunneling). By varying the impact pa- 
rameter and plotting the density contours as a function 
of time, one can easily distinguish between fusion events 
(one fragment in the exit channel) and deep-inelastic re- 
actions (two fragments). 

Let us first consider the case of the neutron-rich system 
^■^^Sn-f-^^Ca. The HF initialization run shows that both 
nuclei are spherical in their ground state, in agreement 
with HFB calculations and with experimental data. At 
energy Ecm. — 130 MeV we find that impact parameters 
< bmax — 4.45 fm result in fusion, while impact param- 
eters b > 4.45 fm lead to deep-inelastic reactions. Using 
the sharp cut-off model, the fusion cross section is given 
by CTfus = T^b'^ax = 62.2 fm^ = 622 mb. This result is dis- 
played in Fig. [31 Similar TDHF runs have been carried 
out at energies E^.m. = 115, 117, 120, 125 MeV, and the 
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corresponding fusion cross sections are also shown in the 
figure. 

When we attempted to calculate fusion cross sections 
for the stable system ^^^Sn+^°Ca we found that the HF 
initialization run correctly predicts a spherical ground 
state for '^'^Ca, but -'^^'^Sn with non-magic neutron num- 
ber = 74 shows a prolate quadrupole deformation 
/32 = 0.054 whereas Skyrme-HFB calculations [s^ 
and experiment show this nucleus to be spherical. The 
unphysical shape is due to the lack of pairing in pure HF. 
Because the shape of the nucleus is very important for 
fusion calculations, we have added BCS/Lipkin-Nogami 
pairing to the HF initialization of ^^^Sn which results 
in an approximately spherical shape. During the sub- 
sequent TDHF run, we have kept the BCS occupation 
numbers frozen. This approximation is necessary, be- 
cause a full 3D-TDHFB code with self-consistent pairing 
for nuclear reactions does not currently exist. 
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FIG. 4. (Color online) TDHF calculation of the total fusion 
cross section at energies above the barrier. At fixed iJc.m. 
energy, the neutron-rich system ^^^Sn+**Ca shows a higher 
cross section than the stable system ^^^Sn+'*''Ca. 

The resulting total fusion cross sections for 
i24Sn-h''°Ca are also shown in Fig. [H The sohd 
black line corresponds to TDHF calculations with added 
pairing (labeled TDHF-hBCS/LN), and the dashed line 
corresponds to pure TDHF which results in a larger 
fusion cross section for the most favorable orientation 
of the (unphysical) deformed ^^*Sn nucleus. We notice 
that if we plot both systems as a function of the 
center-of-mass energy, the more massive neutron-rich 
system shows a larger fusion cross section. Part of 
this is a trivial size-effect which can be taking out by 
appropriate scaling factors. We will comment on this 
issue further when we compare the TDHF cross sections 
to experimental data (Fig. IHlat the end of Section HI). 



B. DC-TDHF heavy-ion interaction potential and 
mass parameter 

In the absence of a true quantum many-body theory of 
barrier tunneling, sub-barrier fusion calculations are re- 
duced to the calculation of a potential barrier between 
the interacting nuclei and a subsequent calculation of 
tunneling through the barrier To date, theoreti- 

cal studies of fusion cross sections are still dominated by 
phenomenological methods such as the coupled-channels 
(CC) approach 37„ 38]. In this approach, one either uses 
empirical ion-ion potentials or one calculates the heavy- 
ion potential with the double-folding method using mea- 
sured nuclear densities for projectile and target [39| . In 
the latter case, one relies on the "frozen density" or "sud- 
den" approximation in which the nuclear densities are 
unchanged during the computation of the ion-ion poten- 
tial as a function of the internuclear distance. The frozen 
density approximation ignores dynamical effects such as 
neck formation during the nuclear overlap. It has been 
demonstrated that for deep sub-barrier energies the in- 
ner part of the potential barrier plays a very important 
role [40,|1. 

While phenomenological methods provide a useful 
starting point for the analysis of multitudinous fusion 
data it is desirable to make contact with the microscopic 
theories of nuclear structure and reactions. We have 
developed a fully microscopic method to extract heavy- 
ion interaction potentials V{R) from the TDHF time- 
evolution of the dinuclear system. In our DC-TDHF ap- 
proach (23| . the time-evolution takes place with no re- 
strictions. At certain times t during the evolution the 
instantaneous TDHF density is used to perform a static 
Hartree-Fock energy minimization. In a typical DC- 
TDHF run, we utilize a few thousand time steps, and 
the density constraint is applied every 20 time steps. In 
this density-constrained method we impose the condition 
that the total proton and neutron density of the dinu- 
clear system at the ion-ion separation distance R{t) is 
equal to the instantaneous TDHF density. We refer to 
the minimized energy as the density constrained energy 
EBc{R{t))- In the DC-TDHF method the ion-ion inter- 
action potential is given by 

ViR) = Euc{R) - Ea, - Ea, , (10) 

where Ea^ and Ea2 are the binding energies of the two 
nuclei obtained with the same effective interaction. For 
the calculation of the ion-ion separation distance R we 
use a hybrid method as described in Ref. ^26*]. At large 
distances where a visible neck allows us to identify two 
fragments we compute it as the distance between the cen- 
ters of mass of the ions. For more compact configura- 
tions, we compute R from the mass quadrupole moment 
Q20 as R = ro \/\Q2o\ where t-q is a scale factor to connect 
the definition smoothly to the large-distance region. 

The interaction potentials calculated with the DC- 
TDHF method incorporate all of the dynamical entrance 
channel effects such as neck formation, particle exchange. 
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internal excitations, and deformation effects. While the 
outer part of the potential barrier is largely determined 
by the entrance channel properties, the inner part of the 
potential barrier is strongly sensitive to dynamical effects 
such as particle transfer and neck formation. 

In Fig. [5] we show the heavy-ion interaction poten- 
tial V{R) for i32Sn-|-48Ca calculated with the DC-TDHF 
method using Eq. (jlOp . The heavy- ion potential depends 
on the center-of-mass energy and is shown at five differ- 
ent Ec.-m. energies. Our results demonstrate that in these 
heavy systems the potential barrier height increases dra- 
matically with increasing energy Ec.m., and the barrier 
peak moves inward towards smaller i?- values. By con- 
trast, DC-TDHF calculations for light ion systems such 
as -'^^0-1-^^0 show almost no energy-dependence even if 
we increase i?c.m. by a factor of four [33|. Even in reac- 
tions between a light and a very heavy nucleus such as 
i6Q_j_208pj-,^ we see only a relatively small energy depen- 
dence of the barrier height and width [26l|. 




R(fm) 

FIG. 5. (Color online) DC-TDHF calculation of the heavy-ion 
potential V(R) for the neutron-rich system ^^^Sn-|-**Ca. The 
energy-dependent potential is shown at five i?c.m. energies. 
The point Coulomb potential is given for comparison. 

In addition to the ion-ion potential, it is also possi- 
ble to compute the corresponding coordinate dependent 
mass parameter M{R) using energy conservation at zero 
impact parameter 

M(«) . , ,n, 

where the collective velocity R is directly obtained from 
the TDHF evolution and the potential V{R) from the 
density constraint calculations. 

In Fig. |6] we show the mass parameter M{R) calcu- 
lated from Eq. pT|) at four different energies -Ec.m.- As 
expected, at large distance R the mass M{R) is equal to 
the reduced mass of the system. At smaller distances, 



when the nuclei overlap, the mass parameter increases 
in all cases. Like the heavy-ion potential, the mass pa- 
rameter shows a strong dependence on energy: at higher 
energies, M{R) rises smoothly with decreasing R. How- 
ever, for the lowest energies which are only a few MeV 
above the corresponding potential barrier, the mass pa- 
rameter shows a pronounced peak. This is due to the 
fact that the velocity R in the denominator of Eq. pT|) 
becomes small as the ions approach each other, reaches 
its minimum at the top of the barrier, and then increases 
again at smaller i?- values. 
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FIG. 6. (Color online) Microscopic mass parameter M{R) 
(in units of the reduced mass p), calculated at several Ec.m. 
energies. 

The total fusion cross section is given by 

^Tfus(^c.m.) = ^(2^ -I- l)Te{Ec.m.) . (12) 

The potential barrier penetrabilities Tg are obtained 
from the numerical solution of the Schrodinger equation 
for the relative coordinate R. Instead of solving the 
Schrodinger equation with coordinate dependent mass 
parameter M{R) for the heavy- ion potential V{R), it is 
numerically advantageous to use the constant reduced 
mass /i and to transfer the coordinate-dependence of the 
mass to a scaled potential U{R) using the scale transfor- 
mation 

dR^i^^y^R. (13) 

Integration of Eq. (fT3|) yields 

R = f{R) ^ R = r\R). (14) 

As a result of this point transformation, the scaled heavy- 
ion potential is given by the expression 

V{R) = V{r\R)) = U{R). (15) 
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In Fig. [7] we compare the original heavy-ion potential 
V{R) (solid lines) corresponding to the mass parameter 
M{R) to the transformed potential U{R) (dashed lines) 
calculated from Eq. ([TSl) which corresponds to the con- 
stant reduced mass /i. We observe that the coordinate- 
dependent mass changes only the interior region of the 
potential barriers, and this change is most pronounced 
at low E'c.in. energies. 
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FIG. 7. (Color online) Solid lines: original heavy-ion po- 
tentials V(-R) corresponding to the mass parameter M{R). 
Dashed lines: transformed potentials U{R) corresponding to 
the reduced mass ^. 

The Schrodinger equation corresponding to the con- 
stant reduced mass /x and the scaled potential U{R) has 
the familiar form 



2fi dR^ 



(16) 

By numerical integration of Eq. using the 

well-established Incoming Wave Boundary Condition 
(IWBC) method 41, 42]wc obtain the barrier penetrabili- 
ties Ti which determine the total fusion cross section (IT^ . 



Fusion cross sections for ^^^Sn+^^Ca and for 
i^*Sn+«Ca 



In Fig. [5] we show the total fusion cross section for 
the neutron-rich system ^■^^Sn-t-'^^Ca. The cross sections 
obtained with unrestricted TDHF runs above the bar- 
rier (see Fig. 2]) are indicated by squares. In addition, 
we display results obtained with the DC-TDHF method 
which can be used at energies Ec.m. below and above the 
potential barriers. The dashed lines correspond to cal- 
culations using the transformed potential U{R) in Fig. [7] 
at three different energies. We observe that at high ener- 
gies, -Ec.m. = 135 — 140 MeV, all of the energy-dependent 



heavy-ion potentials yield the same fusion cross section. 
However, these three potentials yield substantially dif- 
ferent cross sections at i?c.m. energies near and below 
the barrier. In particular, we observe that the heavy- 
ion potential calculated at the lowest possible energy 
Ec.m. = 114 MeV (with a barrier height of 111.9 MeV) 
gives by far the largest fusion cross section at sub-barrier 
energies. According to theory, at a given Ec.m. energy, 
one should use the heavy-ion potential calculated at this 
energy. We call this the "i?c.m. -specific potential". The 
corresponding total fusion cross sections are given by the 
solid black line in Fig. [51 and we can see that these DC- 
TDHF results are in perfect agreement with the unre- 
stricted TDHF runs at above-barrier energies. 




DC-TDHF 

- E - specific potential 

- - potential at E =1 14 MeV 

^ cm. 

- - potential at E =140 MeV - 

^ cm. n 

-- potential at E =180 MeV 



140 



E (MeV) 
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FIG. 8. (Color online) Total fusion cross section for the 
neutron-rich system ^^^Sn+^'^Ca. Compared are DC-TDHF 
calculations using the energy-dependent heavy-ion potential 
in Fig. [7] with unrestricted TDHF calculations (see Fig. [4| at 
energies above the barrier. 

In Fig. [S] we compare our theoretical fusion cross sec- 
tions for i^^Sn-t-^'^Ca and ^^'^Sn-t-^'^Ca to experimental 
data obtained at HRIBF The data for the latter 

system aeree with earlier measurements by Scarlassara 
et al. In order to take out trivial size-effects, the 

HRIBF experimentalists decided to scale their measured 
fusion cross sections and center-of-mass energies as fol- 
lows Defining the size parameter for the composite 
system 



Ro = A 



1/3 



A 



1/3 



the fusion cross section is scaled according to 



(17) 



(18) 



Using the definition of the Coulomb barrier height pa- 
rameter 



(19) 
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the center-of-mass energy is scaled as follows 

Ec.m. ^ Ec.ni./Bo (20) 

This scaling procedure yields the same total fusion cross 
section at high energies for all measured Sn+Ca isotopes. 
Note that when one corrects for nuclear size effects, the 
scaled fusion cross sections for the neutron-rich system 
^^^Sn-f^^Ca are lower than the fusion cross sections for 
the stable ^^^Sn+''*'Ca system, in contrast to the ab- 
solute cross sections plotted in Fig. |4l In the case of 
^^''Sn-|-'*''Ca, we are only able to compare the experi- 
mental data to our unrestricted TDHF runs at above- 
barrier energies (see Fig. H]). As we pointed out in Sec- 
tion III. A, these runs were initialized with BCS/Lipkin- 
Nogami pairing for ^^^Sn, with the BCS occupation num- 
bers kept frozen during the collision. Unfortunately, 
this approximation cannot be utilized in the DC-TDHF 
method because the static HF solution coupled with a 
constraint on the instantaneous TDHF density for the 
combined system requires the reevaluation of the occu- 
pation numbers for the lowest energy solution. The good 
agreement between the measured fusion cross sections 
and our theoretical results is quite remarkable because 
the only input in TDHF and DC-TDHF is the Skyrme 
N-N interaction with no adjustable parameters. 




FIG. 9. (Color online) Comparison between scaled experimen- 
tal fusion cross sections [31 and theoretical results obtained 
with the DC-TDHF method and with unrestricted TDHF. 



IV. SUMMARY 

In this paper, we present dynamic microscopic calcula- 
tions of pre-equilibrium giant dipole resonance excitation 
and fusion in the neutron-rich system ^^^Sn-|-''*Ca at en- 
ergies near the Coulomb barrier, and we compare photon 
yields and total fusion cross sections to those of the stable 
system ^^''Sn-|-'"'Ca. The numerical calculations are car- 
ried out on a 3D lattice using both the Time-Dependent 
Hartree-Fock (TDHF) method and the recently devel- 
oped Density Constrained TDHF method (DC-TDHF). 

Regarding pre-equilibrium GDR excitation, we demon- 
strate that the peak of the dipole radiation spectrum 
occurs at a substantially lower energy than expected for 
an equilibrated system, thus reflecting the very large pro- 
late elongation of the di-nuclear complex during the early 
stages of fusion. 

In addition, we study total fusion cross sections for 
both systems. First, we carry out unrestricted TDHF 
runs at energies above the barrier (no tunneling). Sec- 
ondly, we calculate the heavy-ion interaction poten- 
tial V{R) and the coordinate-dependent mass parame- 
ter M{R) utihzing the DC-TDHF method. Our resuhs 
demonstrate that in these heavy systems the potential 
barrier height increases dramatically with increasing en- 
ergy Ec.m., and the barrier peak moves inward towards 
smaller i?- values, in contrast to light ion systems such as 
16q_|_16q .^yjjich show almost no energy- dependence. 

By numerical integration of the Schrodinger equation 
for the relative coordinate R with the Incoming Wave 
Boundary Condition (IWBC) method, we calculate to- 
tal fusion cross sections both below and above the bar- 
rier. We emphasize that the only input in our micro- 
scopic calculations is the Skyrme N-N interaction, with 
no adjustable parameters whatsoever. We find that our 
theoretical fusion cross sections for both systems agree 
reasonably well with recent data measured at HRIBF. 
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